Load libraries

library(tidyverse)
library(plotly)
library(readxl)

Set WD and read Excel file

setwd("~/GitHub/DSCI605/Mod9")
land = readxl::read_xlsx('Crop_Range_GOES0901_CountJday.xlsx')

Create plot

#Histogram
plot_ly(land, x = ~FDCount,nbinsx = 10) %>%
  filter(jday %in% 233) %>%
  add_histogram(marker = list(color = "teal",
                              line = list(color = "darkgray",
                                          width = 2)),
                name = "Julian Day = 233") %>%
  layout(title="FDCount Histogram at Julian Day 233",
         yaxis=list(title = "Count"),
         xaxis=list(title = "FDCount")) %>%
  layout(showlegend = TRUE)

Another way

plot_ly(land, x = ~FDCount,nbinsx = 10, type ="histogram") %>%
  filter(jday %in% 233) %>%
  layout(title="FDCount Histogram at Julian Day 233",
         yaxis=list(title = "Count"),
         xaxis=list(title = "FDCount")) %>% 
  layout(showlegend = TRUE)